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H ' Abstract 

We present CO interferometer observations of the barred galaxy NGC 4303 (M61). This galaxy has 
a strong gas eoncentration at the central region and prominent offset ridges at the leading sides of the 
bar. Sharp velocity gradients are apparent across the ridges. The brightness temperature in the ridges 
is low, indicating the existence of unresolved molecular clouds. Additionally, an analysis of the newborn 
stellar clusters revealed in HST images suggests that the parent molecular clouds had masses of IQ'^^^Mq. 
The observed shear velocity gradient across the ridges is too small to break up giant molecular clouds. 
' Therefore, the clouds are likely to survive passage through the ridges. We discuss a cloud orbit model 

\ in a bar potential for understanding the gas distribution, dynamics, and star formation in NGC 4303. 

^SJ ■ The model reproduces the narrow offset ridges and sharp velocity gradients across the ridges, although 

I no shock is associated with the ridges. We discuss cloud-cloud collisions (and close interactions) as a 

, possible triggering mechanism for star formation. The newborn stellar clusters in NGC 4303 are located 

predominantly at the leading sides of the offset ridges. Cloud orbits are densely populated in the region 
where the clusters are found, suggesting a high coUisional frequency and possibly a high rate of triggered 
' star formation. Cloud-based dynamics is less dissipative than smooth hydrodynamic models, possibly 

extending the timescales of gas dynamical evolution and gas fueling to central regions in barred galaxies. 
Q i Key words: galaxies: individual(NGC 4303) - galaxies: ISM - galaxies: kinematics and dynamics 

, 1. Introduction associated gas kinematics in bars. 

The large-scale distribution of star-forming regions and 

- ^ , Despite numerous observational and theoretical studies, young stars provides a clue to understanding gas dynam- 

■ the nature of interstellar medium (ISM) gas dynamics re- ics. Star formation triggered by hydrodynamic shocks 

] mains uncertain (Combes 1996). The dynamical response would be localized in narrow shocked regions. Cloud-cloud 

of the gas depends significantly on whether the ISM is a interactions could potentially occur everywhere in galactic 

continuous hydrodynamic fluid or a set of discrete molecu- disks and cause star formation. We compare CO observa- 

lar clouds. The hydrodynamic fluid produces shocks along tions of the barred galaxy NGC 4303 to the cloud-based 

spiral arms or bars in galaxies, while the molecular clouds gas dynamics model. The cloud-based model naturally 

could ballistically pass the spiral/bar structures without explains the distribution of young stellar clusters, 

encountering a shock. Such differences must change not Figure 1 shows the B and Ha band images of NGC4303; 

only the dynamical response of the gas (e.g. shock, dissi- most HII regions are associated with outer spiral arms, 

pation), but also the evolution timescale (e.g. gas inflow) but almost no apparent Ha emission is found in the bar. 

and star formation. Similar trends have been found in many other barred 

The Galactic ISM is concentrated in discrete molecu- galaxies (Downes et al. 1996; Regan et al. 1996; Garcia- 

lar clouds (Scoville, Sanders 1987), and supports the dy- Barreto et al. 1996; Sheth et al. 2000; 2002). The ab- 

namics of slightly coUisional clouds. Individual molecular sence of active star formation in bars is often attributed to 

clouds are rarely resolved in external galaxies; however, strong shear velocity in the bar (Reynaud, Downes 1998). 

there is a general belief that only continuous hydrody- However, we point out below that the shear in NGC 4303 

namic models could reproduce the narrow gas/dust lanes is too small to break up dense molecular clouds, 

observed in barred galaxies, since hydrodynamic fluids NGC 4303 hosts remarkable central activity: a nuclear 

easily produce narrow shock lines along bars. We demon- secondary bar (Colina, Wada 2000), prominent UV emis- 

strate that an alternative scenario, cloud-based gas dy- sion around the secondary bar (Colina et al. 1997), and 

namics, also reproduces the narrow gas distribution and an active galactic nucleus (AGN; Colina, Arribas 1999). 
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It has been a prime example of the starburst-AGN con- 
nection (Colina et al. 1997; Colina, Arribas 1999; Colina, 
Wada 2000; Jimencz-Bailon ct al. 2003). Schinncrer et 
al. (2002) discussed gas fueling toward the nucleus due 
to the secondary bar based on their CO data. Our fo- 
cus, however, is on the gas dynamics and star formation 
around the offset ridges in the primary bar. 

We describe the observations and data in §2, and re- 
sults in §3. We introduce a cloud-orbit model in the bar 
and show its ability to reproduce the observed CO dis- 
tribution and velocity structure (§4). We discuss star 
formation along the bar using the cloud-based model 
(§5). Conclusions are summarized in §7. This paper 
is the sixth in a series on the Virgo CO survey at the 
Nobeyama MiUimeter Array (Sofue et al. 2003a; 2003b; 
2003c; Onodera et al. 2004; Nakanishi et al. 2005). 

Figure 1 

Table 1 

2. Observational Data 

2.1. NMA CO (1-0) Observations 

We made aperture synthesis observations of the Virgo 
cluster galaxy NGC 4303 in the CO J = 1 — line emis- 
sion using the Nobeyama Millimeter Array (NMA). The 
observations were made between 1999 December and 2000 
February with a single pointing center at (aigso, <5i95o) 
= (12'^19'^2P.6, -h4°45'3".0). Six 10 m telescopes pro- 
vide the FWHP of about 65" at 115 GHz. We used 
three available array configurations (AB, C, and D con- 
figurations). The typical system noise temperature was 
about 400 K in single sideband. Digital spectro-correlators 
(Okumura et al. 2000) have two spectroscopic modes; we 
used the mode covering 512 MHz (1331 kms~^) with the 
2 MHz (5.2kms~^) resolution. We observed the quasar 
3C273 every 20 minutes for gain and bandpass calibra- 
tions. The absolute flux of 3C273 was measured against 
Uranus and/or available flux-known quasars; it was 14.2Jy 
for 1999 December, and 11.4Jy for 2000 January and 
February. The uncertainties in the flux scale are about 
20%. 

Raw visibility data were calibrated for pass- 
band response and complex gain variation with the 
NRO/UVPROC-II package (Tsutumi et al. 1997). The 
data were binned to 4 MHz resolution. Each channel 
data was deconvolved with the CLEAN procedure using 
the NRAO/AIPS package. We applied both natural and 
uniform weightings for the visibility data and obtained 
2". 4- and 1". 3- resolution maps, respectively. We refer the 
2".4-resolution map as "map". The parameters of the 
observation and data cubes are given in Table 2 and 3, 
respectively. 

Table 2 

Table 3 

Figures 2 shows zeroth- and first-moment maps, and 
synthesized beam map. Figure 3 and 4 show the velocity 



channel maps and position-velocity (PV) diagrams, re- 
spectively. Primary beam attenuation is not corrected 
for these maps. In the synthesized beam map (Figure 
2), prominent sidelobcs run in the north-south direction 
and are aligned with the observed emission distribution. 
These sidelobes made the deconvolution less obvious; we 
arranged restricted CLEAN boxes and checked that the 
2nd- and 3rd-order sidelobes are also reduced successfully 
in the deconvolution. We estimated the recovered flux by 
comparing our data with a single-dish observation at the 
FCRAO 14 m telescope (Kenney, Young 1988). Our map 
recovers 97% of the total flux in the central 45" (with the 
uncertainties of 20%) . 

Figure 2 

Figure 3 

Figure 4 

2.2. Supplied Optical/Near-IR Data 

We obtained optical and near-infrared images from 
the archives and literature through the NASA/IPAC 
Extragalactic Database (NED): the B-band image from 
the Digitized Sky Survey, Ha-image (Figure 1 and 5) from 
Koopmann, Kenney, and Young (2001), and A'-band im- 
age (Figure 5) from MocUenhoff and Heidt (2001). We 
also obtained HST images in the F450W and F814W- 
filters (P.I. Smartt) through the HST archive at the 
Canadian Astronomy Data Center (CADC). The F450W 
and F814W images are shown in Figure 6. 

The astrometry of the B- and Ha-band images was 
determined using the USNO-A2.0 catalog (Zacharias et 
al. 2000). The ivT-band and HST images do not contain 
enough cataloged stars; we shift the central compact emis- 
sion peak in each image to the dynamical center (§3.3). 
The centroids of compact sources across these images co- 
incide with those in the lower-resolution _B-band image. 
The uncertainty of the absolute position in each image is 
about 0".5. 

Figure 5 

Figure 6 

3. Results 

Figure 2 shows two remarkable features in the gas distri- 
bution: (1) a central concentration within r ~ 10" (780pc, 
hereafter the central disk), and (2) offset ridges extending 
from the central concentration out to r ~ 30" (2.3kpc). 
The gas is rotating clockwise, assuming trailing spiral 
arms. Below, we obtain the parameters of the gas disk 
and discuss the two features. 

3.1. Global CO Properties 

Figure 7 shows the radial profile of the CO inte- 
grated intensity Ico^V cos(i) and cumulative flux ^co- 
Corrections for the primary beam response and inclination 
were applied. We adopted the position angle P. A. = —45°, 
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inclination i = 30° (see §3.2), and the CO-to-H2 con- 
version factor Xco = 1-8 x 10^°cm^^ [Kkms^^]^^ (Dame 
et al. 2001). The surface density is ahnost constant at 
~ 5OOM0pc~^ from the center to r 4" (310 pc), and 
decreases exponentially with a scale length of ~ 4.5" 
(350pc). The total detected gas mass is 8.4 x IO^Mq 
within the radius 13" (Ikpc), and is 1.0 x IO^AIq within 
30" (2.3kpc). 

Figure 7 

3.2. Central Concentration 

The CO emission is very concentrated in the central disk 
(r < 10") in Figure 2. The velocity field is the typical "spi- 
der diagram" . The gas is following almost pure circular ro- 
tation. Spiral arms are curling along the outer edge of the 
central disk (Figure 8). They are associated with slight 
perturbations in the velocity field at the outer edges of 
the central disk. We obtained kinematic parameters from 
the velocity-field assuming pure circular rotation. The dy- 
namical center, position angle, inclination, and recession 
velocity were determined using the AIPS/GAL package 
(Table 4). Since the perturbations due to the spiral arms 
are small, the errors due to the perturbations should be 
small. The derived parameters are consistent with those 
from previous studies (Colina, Arribas 1999; Schinnerer et 
al. 2002). 

Figure 8 

3.3. Offset Ridges 

Two gas ridges run out from the central disk in the 
north and south directions (Figure 2). Figure 5 (left) com- 
pares the gas distribution (grayscale; CO image) with the 
stellar bar (contours; K-hand image). The two ridges are 
at the leading sides of the stellar bar. This type of ridge is 
referred to as "offset ridges" , and is often found in barred 
galaxies (Ishizuki et al. 1990; Sakamoto et al. 1999; Sheth 
et al. 2002). The gas ridges coincide with dust lanes in 
Figure 6. There are gas concentrations at the outermost 
ends of the detected CO ridges (Figure 2), which coin- 
cide with the stellar bar ends. Such gas concentration is 
common in barred galaxies. Note that no primary beam 
correction has been applied to this map; the fluxes of these 
gas concentrations are therefore about 2-timcs larger than 
the contour values shown in the figure. 

Sharp velocity gradients exist across the offset ridges; 
in the channel maps (Figure 3), the southern CO ridge 
shifts east (toward left) as the velocity decreases, and the 
northern ridge shifts west (right) as the velocity increases. 
The isovelocity contours run along the offset ridges in the 
velocity field map (Figure 2). The velocity differences 
across the ridges, ^ 4", is about ^ 20kms^^ (i.e. three 
channels) in projection on the sky. 

The typical brightness temperature is Tb 1 K along the 
offset ridges (Figure 3). Assuming an excitation tempera- 
ture Tex = lOK for optically-thick molecular gas (Scoville, 
Sanders 1987) and the cosmic microwave background tem- 
perature of TcMB = 2.7K, the area filling factor within 
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the beam (2.4" - 200 pc) is / ~ 0.1 . Thus, the emis- 
sion arises from unresolved clumpy structures, presumably 
molecular clouds. The average gas surface density in the 
beam is typically Sgas — 15OM0pc~^. Therefore, molec- 
ular clouds are dense with an average surface density of 
Smc ~ Sgas// ~ IOOOMqPC"^. This is much denser than 
17OM0pc~^ for clouds in the Galactic disk (Solomon et 
al. 1987), but less dense than 25OOM0pc~^ for clouds in 
the Galactic center (Oka et al. 2001). Note that a veloc- 
ity dilution might underestimate the value / by a factor 
of 2, although it does not affect the following discussion. 
The typical velocity width of the offset ridges is about 
20kms^^ (Figure 3), and that of molecular clouds in the 
Galaxy is 9kms-i (FWHM; Scoville, Sanders 1987). 

3.4. Star Formation around Offset Ridges 

Evidence of star formation (SF) around the offset ridges 
is inferred from many bright points in the HST images 
(Figure 6), although active SF is not obvious in the Ha 
image (Figure 5 right). Each of the bright points is very 
likely a stellar cluster. Their luminosities indicate that the 
masses are 1O^~^M0. Assuming an extinction Ay = Imag 
and a color excess 0.71 (Cardelli et al. 1989), the typical 
intrinsic color in F450W-F814W is as blue as —0.7 mag 
(AB magnitude), indicating their young ages, ^ 10^ yr 
(see Appendix). Because of the large scatter in color 
(±0.5 mag), we could not detect any significant system- 
atic color change as a function of the distance from the 
ridges. 

The clusters are distributed predominantly at the lead- 
ing side of the gas ridges, i.e. at the downstream side as- 
suming that the corotation radius is around the bar end. 
No obvious cluster is found at the opposite (upstream) 
side, which is especially notable for the southern ridge. 
Therefore, the star formation was initiated during, or af- 
ter, passage of the parent molecular cloud across the offset 
ridges. Some clusters are found far from the gas ridges. 
The rotation timescale around the galaxy, lO'^yr, is ten- 
times longer than the cluster age, ^ 10^ yr. Thus, those 
clusters should have been formed not on the ridges, but 
long after passage of the cloud across the ridges. 

The conversion of molecular gas to stars in a molecular 
cloud is normally a very inefficient process in the Galaxy 
with a conversion efficiency of 0.01-0.1 (Evans, Lada 1991; 
Lada 1992). Consequently, the parent molecular clouds of 
the stellar clusters should have had masses of 1O'*~^M0, 
which arc typical for Galactic molecular clouds (Scoville, 
Sanders 1987). 

3.5. Comparison with Other CO Maps 

Three CO maps of NGC 4303 have been published 
from mm- wave interferometer observations, at the Owens 
Valley Radio Observatory (OVRO; Schinnerer et al. 2002), 
at the Nobeyama Millimeter Array (NMA; Sofue et 
al. 2003a, and this paper), and at the Bcrkeley-IUinois- 
Maryland Association (BIMA; Heifer et al. 2003). 

The BIMA map has three-times lower sensitivity of 
~ 47mJy and resolution of ~ 6", but includes zero-spacing 
baselines (extended components) from single dish obser- 
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vations. It is thus good for tracing the large-scale distri- 
bution of molecular gas. The BIMA map shows significant 
CO emission along the outer spiral arms. The brightness 
is comparable between the offset ridges and spiral arms, 
though star formation is more active in the arms (Figure 

The OVRO map has similar sensitivity (~ 20mJy) and 
resolution (~ 2") to our map. The CO distribution and 
velocity structure are consistent between the OVRO and 
NMA maps. The recovered flux compared with the single- 
dish data (Kenney, Young 1988), however, is 2-times 
higher in NMA than in OVRO. We compared the two 
data cubes^, and found that the NMA flux is ~ 1.4 times 
larger in spectrum and ^ 2 times greater in integrated 
intensity. The lower OVRO flux probably results from a 
combination of the assumed calibrator flux (~ 10 Jy vs 

14.5 Jy) and a lack of short-spacing uv coverage in the 
OVRO data. This difference does not affect the discus- 
sions made in this paper and in Schinnerer et al. (2002). 

4. Gas Dynamics in NGC 4303 

We use a cloud-orbit model to interpret the gas dy- 
namics in NGC 4303. To derive the stellar orbits in the 
bar, the equations of motion for a test particle are solved 
in a bar potential (Contopoulos, Papayannopoulos 1980; 
Binney, Tremaine 1987). Wada (1994) included the 
damping force term, "— 2Ai?", in the equations of mo- 
tion, the force negatively proportional to the velocity R, 
and obtained the orbits of coUisional gas clouds. This 
model has provided successful interpretations for observed 
barred and spiral galaxies (Sakamoto et al. 2000; Koda et 
al. 2002; Onodera et al. 2004). Figure 9, 10, and 11 show 
the model gas orbits, density and velocity field maps, and 
position- velocity diagrams. The gas and bar are rotating 
clockwise. The model is scalable; it fits to NGC 4303 when 
the units are 300pc in length and 160kms~^ in velocity. 
The pattern speed and bar radius (corotation radius) are 
scaled to 53kms~^ kpc~^ and 3kpc, respectively, which is 
consistent with the bar radius in Figure 1. 

Figure 9 

4-1. Offset Ridges and Central Spirals 

Figure 10 shows the projection of the gas orbits with a 
position angle of —45° and an inclination of 30° (Table 4). 
The stellar bar runs vertically. The density map (middle) 
reproduces the observed offset ridges and the central spi- 
ral arms in Figure 2. To calculate the density, we used the 
timescale of passage along the orbit as a weighting func- 
tion. The offset ridges and central spiral arms appear to 
be dense for two reasons: (1) the gas orbits are crowded at 
the leading side of the bar, and (2) passing the apocenter 
of an orbit takes longer time than passing the pericenter. 

There is a common expectation that the offset ridges are 
produced by hydrodynamic shock; the tangential orbital 
velocity of the gas is damped by the shock, and the gas 

^ OVRO data cube was kindly provided by Eva Sciiinnerer. 



falls toward the galactic center along the offset ridges. The 
orbit model, however, indicates that no shock is required 
to form the dense and narrow offset ridges. In fact, the 
sharp turns are apparent in the rotating frame with the 
bar (Figure 9 left), but disappear in the rest frame (right). 
The tangential velocity becomes very close to the rotation 
speed of the bar on the offset ridges, which produces the 
apparent turns. This fact is particularly important when 
the orbits of young stars born recently in molecular clouds 
are discussed. We note that this argument does not ex- 
clude the possibility of shock as a secondary effect. 

Figure 10 

The velocity field (Figure 10 right) is consistent with 
the observations (Figure 2 middle). The major axes of 
elongated gas orbits are aligned with that of the galaxy at 
the central disk, and therefore, the isovelocity contours are 
similar to those of pure circular rotation (spider diagram). 
Slight deviations from pure circular rotation are apparent 
along the spiral arms on the central disk; we found simi- 
lar deviations in the observed velocity field (Figure 2 mid- 
dle). At the outer part, the isovelocity contours nm along 
the offset ridges as found in the observations. The con- 
tours parallel to the offset ridges indicate abrupt velocity 
changes across the ridges, although no shock is associated. 

Figure 11 

Figure 11 shows position- velocity (PV) diagrams with 
three different position angles (PA). The diagram with 
PA = —45° corresponds to NGC 4303 and is similar to 
the observed one. The model reproduced the observations 
in space (maps) and velocity (PV diagram), indicating 
that the orbit model represents the gas motions in NGC 
4303. The other two diagrams, i.e. PA = 0° and -90°, 
are projections from the edge-on and pole-on directions 
of the bar, respectively. The same gas distribution shows 
quite different patterns; our determination of the PA and 
inclination in §3.2 were suitable for NGC 4303. 

4.2. Cloud-cloud Collisions as a Damping Force 

The damping force term "— 2Ai?" provides the damp- 
ing timescale, tdamp = 1/A, in which random motions 
decay. For the above model, we defined the damp- 
ing coefficient A = A/kq = 0.1 as in previous studies 
(Wada 1994; Sakamoto et al. 2000). The orbits are insensi- 
tive to the change of A by a factor of a few (Wada 1994). 
We discuss its physical cause in an order-of-magnitude 
approximation. Assuming a constant rotation velocity 
V = 160kms-i at a radius of Ikpc as for NGC 4303, the 
epicyclic frequency becomes kq ~ 200kms~^ kpc~^. The 
damping timescale is then tdamp = 1/Ako ~ 5 x 10^ yr for 
A = 0.1. 

We have shown the presence of dense molecular clouds 
in NGC 4303 (§3.3). The damping timescale is close 
to the collisional timescale of clouds calculated as fol- 
lows. NGC 4303 has a gas mass of IO^Mq between 
radii of 1 — 2kpc (Figure 7). Assuming a cloud mass 
Mmc = IO^Mq and diameter D = 20pc, the surface num- 
ber density of molecular clouds is iVMc = 100 kpc^^. If 
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the cloud-cloud velocity dispersion is cr = lOkms^^ and 
the galactic molecular disk is thin (2-D) in the central re- 
gion, the coUisional timescale between molecular clouds 
becomes tcoi = I/NmcDct ^ 4 x 10 ""yr. Close encounters 
would occur on a slightly shorter timescale but on the or- 
der of lO*" yr. These are close to the damping timescale 
for A ~ 0.1. Cloud collisions could be the cause of the 
damping force. 

Ram pressure from ambient gas is negligible. Using the 
density of the ambient gas p, the relative velocity between 
a cloud and the ambient gas v, the cloud cross section 
A, and the mass AImc, the acceleration on the cloud is 
ftram = Pv'^A/Mmc- The velocity difference would decay 
in ^ram = ''^'/'^ram- If the ambient gas is atomic hydrogen 
with a density of Icm"'^ and a velocity of v = lOkms"^, 
iram ~ 10^ yr, which is much longer than icoi- 

4-3. Shear Velocity and the Survival of Molecular Clouds 

Downes et al. (1996) discussed that molecular clouds 
might be destroyed by strong shear around offset ridges. 
Assuming that the shapes and orientations of the gas or- 
bits in Figure 10 left are all correct, molecular clouds have 
only the velocity components tangential to our line-of- 
sight before the ridges, and only the line-of-sight veloc- 
ity components after the ridges. The velocity difference 
across the ridges is then ^ 40kms~^ from the observed 
velocity difference of 20kms~^ and an inclination of 30°. 
The width of the ridge is ~ 300 pc (~ 4"), and thus the 
velocity gradient (shear) is about 0.13kms~^ pc~^. 

The surface density of a molecular cloud in NGC 4303 
is about Smc = IOOOMq pc~^ (§3.3). Assuming a cloud 
radius R = 10 pc, the escape velocity from the cloud, 
Vesc = i^nGTjMcRy^^ , is 16kms~^. Thus, a sharp ve- 
locity gradient of 1.6kms~^ pc^^ is necessary to destroy 
the molecular cloud. This is an order of magnitude above 
the observed shear velocity gradient in NGC 4303. Even 
for some molecular clouds with lower surface densities of 
lOOA'/0 pc"^ (typical for the Galactic molecular clouds), 
the escape velocity is still higher than the observed shear 
velocity. Thus, shear is unlikely to destroy molecular 
clouds. The clouds are likely to survive passage across 
the ridges. 

5. Collision-Induced Star Formation in the Bar 

The existence of molecular clouds is necessary but 
not sufficient for star formation (Mooney, Solomon 1988; 
Scoville, Good 1989; Fukui et al. 1999). There is con- 
siderable evidence that cloud-cloud collisions may initiate 
star formation (Loren 1976; Scoville et al. 1986; Odenwald 
et al. 1992; Hasegawa et al. 1994). The efficiency of col- 
lisional dissipation (§4.2) also supports the cloud-based 
mechanism for triggering star formation. We explore the 
cloud-cloud collision (or close encounter) model to account 
for star formation. 

5.1. Star Formation at the Leading Side of Ridges 

The damping timescale of random motions is close to 
the rotation timescale of 10^ yr around the galactic cen- 



ter (§4.2). Thus, the dissipation is not very effective, and 
newborn stars should follow similar orbits to those of the 
parent molecular clouds (Figure 9) within about one rota- 
tion timescale. Comparing the distribution of the young 
clusters (Figure 6) to the orbits (Figure 10 left), we find 
that the majority of the clusters are aligned on the out- 
ermost orbits in Figure 10. Hence, most star formation 
occurs at the outer ends of the offset ridges. The same 
consideration can be applied to other barred galaxies (see 
Sheth et al. 2002); the outer ends of offset ridges are the 
primary loci of star formation in the bars. 

We discussed in §3.4 that some of the stellar clus- 
ters have large offsets from the gas ridges (Figure 6). 
Considering the isochrones along the orbits (Figure 12) 
and the young cluster ages ^ lO^yr, they should have been 
born long (10'' yr) after the parent clouds cross the ridges. 
Star formation is, therefore, triggered even after (but not 
before) the offset ridges. This cannot be explained by the 
star formation induced by galactic shocks. In Figure 9 
left, the separation between orbits is remarkably narrow 
at the leading side of the ridges. The number density of 
clouds would therefore be larger at the leading side, and 
possibly result in more frequent cloud collisions (or close 
encounters). The collision model explains the star forma- 
tion after the ridges. 

Figure 12 

The presence of HIT regions at the leading side of 
the ridges has been a problem for hydrodynamic models 
(Sheth et al. 2002). They predict that the gas looses the 
tangential velocity at offset shocks, and falls toward the 
galactic center along the offset ridges (Regan et al. 1999, 
see their Figure 9). If star formation occurs due to the 
shocks, i.e. after the tangential velocity is damped, the 
newborn stars should also fall along the offset ridges and 
cannot exist at their leading side (Sheth et al. 2002). The 
orbit model, however, naturally explains the co-motion of 
molecular clouds and stars, and the star-forming regions 
at the leading side. 

Young stars are bright in UV emission for 10^ yr after 
birth (Iglesias-Paramo et al 2004), which is close to the 
rotation timescale. If the young stars follow the gas or- 
bits for a single rotation timescale and if the outer edges 
of the ridges are the primary loci of star formation, the 
outermost orbits in Figure 10 left would stand out in UV 
emission from young stars. This is the case for NGC 4303 
in the GALEX/UV image (Gil de Paz et al. 2006 in prepa- 
ration). 

6. Cloud Model vs Hydrodynamic Model 

Most hydrodynamic models assume an isothermal gas 
with a sound speed of ~ lOkms"^, which is similar to 
the observed cloud-cloud velocity dispersions. The corre- 
sponding effective temperature is about 10^ K. far higher 
than the typical kinetic temperature of molecular gas 
(~ 10 K). Therefore, the hydrodynamic models represent 
not the diffuse molecular gas. but an ensemble of a signif- 
icant number of molecular clouds. The difference between 
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the two models arises from the number of clouds being 
considered. Cloud-based dynamics assumes fewer clouds, 
and thus, would be less dissipative than smooth hydro- 
dynamic models. The efficiency of gas fueling to galactic 
centers could be lower than what we expect from hydro- 
dynamic simulations. 

Obviously, there is diffuse, presumably continuous, gas 
around the discrete clouds. The gas is cycling between 
the two phases due to cloud formation, stellar feedback, 
ram pressure stripping of the cloud envelope, etc. This 
cycling would affect the gas dynamics. To discuss these 
effects, molecular clouds have to be resolved in galaxies 
with higher resolution observations (e.g. ALMA) and sim- 
ulations (e.g. Wada, Koda 2001). 

7. Conclusions 

We observed the central 5 kpc of the barred galaxy NGC 
4303 in the CO (1-0) line with the Nobeyama Millimeter 
Array. We discussed the gas distribution and dynamics, 
and compared them with young stellar cluster distribution 
in HST images. The main conclusions are: 

1. NGC 4303 has a central gas concentration (disk) and 
offset ridges along its stellar bar. The central disk follows 
pure circular rotation with slight perturbations due to spi- 
ral arms on the disk. Sharp velocity gradients arc observed 
across the offset ridges. 

2. The brightness temperature in the 2". 4- resolution 
map is much lower than the typical excitation temperature 
of molecular gas in the Galaxy. Assuming that the gas is 
optically thick in CO (1-0), the beam-filling factor is about 
0.1, indicating the existence of dense molecular clouds (~ 
IOOOMq pc"2) in the offset ridges. 

3. The velocity gradient across the offset ridges is about 
0.13kms"i pc"^ This is an order of magnitude smaller 
than the gradient necessary to break up existing molecular 
clouds. Therefore, the clouds survive after passing the 
offset ridges. 

4. Many young stellar clusters are found in HST images ^ 
although star formation is not apparent in the Ha image. 
The clusters have ages of ^ 10^ yr and masses 1O^~'*M0. 
If the star-formation efficiency is 0.01-0.1, their parent 
molecular clouds should have had masses of lO^^^Af0. 

5. The young clusters are distributed predominantly at 
the downstream side of the offset ridges. Some of them 
are located far from the offset ridges and, therefore, are 
born long after their parent clouds pass the ridges. 

6. We used a cloud orbit model to simulate the gas 
dynamics. The dynamics due to slightly-coUisional clouds 
reproduces the central gas disk, narrow offset ridges, and 
their associated velocity structure. The offset ridges are 
not necessarily associated with shock. 

7. The strength of viscous (damping) force is a param- 
eter of choice in the cloud orbit model. The best-fit value 
gives a timescale of 5 x 10^ yr for the decay of random 
motions, which is close to the timescale of cloud-cloud 
collisions or close interactions in NGC 4303. 

8. Cloud collisions naturally explain the formation of 
the young clusters far from the offset ridges. Cloud-orbits 



are densely populated at the leading side of the ridges, 
indicating a high collision frequency and star formation. 

Appendix 1. Color Evolution of a Stellar Cluster 

Figure 13 shows the model color evolution of a stel- 
lar cluster. Three photometric synthesis models are com- 
pared: Bruzual and Chariot (1993, hereafter, BC93), Fioc 
and Rocca-Volmerange (1999, Pegase), and Leitherer et 
al. (1999, Starburst99). We used the Salpeter initial mass 
function (IMF) with a mass range 0.1 - IOOMq for BC93, 
the Scalo98 IMF with a range 0.1 - 120^© for Pegase, and 
the IMF having indices of 1.3 for 0.1 - O.SM© and 2.3 for 
0.5 - IOOMq for Starburst99. The F450W-F814W color 
drops sharply around 7 x 10^ yr for all three models. The 
uncertainties are a factor of 2 in time and 0.1 mag in mag- 
nitude. The stellar clusters bluer than -1-0.0 mag (§3.3) 
are, therefore, younger than ^ 10^ yr (and 2 x 10^ yr at 
the oldest). 

Figure 13 
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Fig. 1. Le/t; i?-band image of NGC 4303 from the Digitalized Sky Survey. Right: Ha image from 
Koopmann, Kenney, and Young (2001). The circles represent the primary beam size of the ^^CO(J =1 — 0) 
observations (65" HPBW). 




UZ]4tJ! UM SIA £4,1 UJb 




tI£1»jD 55,1} 5*^ ».« 

RioHT*scEr(siori ijsoonii 



■IE -1Q -fl D ft 1<t 19 
AflCSEC 



GrEyM:nlrHlLU!ranjiB.[>.ao lii.ai jybeam'Jtinfi i£mBcaef!lJ!lfanoE-lilT4.S ttims km/k L*1(S = tflMEtJl " (-1 -1, 1, 3, S, 5, 5,7. 

uvt=i53i "lii.i.a. ii.i. e.T. a, ». i») i.*v* = i;i-iiD. issS. iwo. ism, i"" — ~ ~ 



.^5»B.•lfm.^^m) e, i, loi 



Fig. 2. CO(l-O) zcroth- and first-moment maps of NGC 4303, and synthesized beam map. A cUp at the 
2a level was used to make these maps. The circles represent the size of the primary beam. No primary 
beam correction was applied. 
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Fig. 3. CO(l-O) velocity ciiannel maps of NGC 4303 in the same region as in Figure 2. The velocities 
at the channels arc shown at the lower right corners. No primary beam correction was applied. The 
ImJybeam^^ corresponds to 0.017K in brightness temperature. 
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Fig. 4. Position-velocity diagrams with PA = —44.45 and slit widths of 30" {upper panel) and 3" {lower 
panel). The contours arc 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0 times the peak intensity. 



10 




i:iTI7 » 5& 5H S iZ 12 2T5T M 59 54 53 



Fig. 5. K-hand and Ha images of tlie central 90" x 90" region of NGC 4303. Upper-left: A'-band image. 
Lower-left: if-band contours on CO grayscale image. Contours are 0.008, 0.012, 0.018, 0.027, 0.041, 
0.061, 0.091, 0.14, 0.21, 0.31, 0.46, 0.69, 1.0 times the peak if-band intensity. Upper-right: Ha image. 
Lower-right: Ha contours on CO grayscale image. Contours are 0.031, 0.063, 0.125, 0.25, 0.5, 1.0 times 
the peak Ha intensity. The circles represent the field-of-view (FWHM) of CO observations. 



11 




Fig. 6. HST images of tlic central 90" x 90" region of NGC 4303, the same region covered in Figure 1. 
Upper panels: HST F450W band images. Lower panels: HST F814W band images. Right panels are the 
central vertical strips of the left images, and CO contours arc overlaid. Contours arc 0.1, 0.2, 0.3, 0.4, 
0.5, 0.6, 0.7, 0.8, 0.9, 1.0 times the peak CO intensity. The circles represent the FoV (FWHM) of CO 
observation. The bright points across the images are very likely stellar clusters. They are distributed 
predominantly at the leading (downstream) side o^^e gas ridges, i.e. the eastern side of the southern 
ridge and the western side of the northern ridge. 
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Fig. 7. Radial profiles of tlie integrated- intensity (inclination corrected). The physical scales at 
the top and right axes were calculated using the distance of 16.1 Mpc and a conversion factor of 
-^co = 1-8 X 10^" cm~^[Kkms~^]~^. We assumed Mgas ~ lAlM^^, taking He and the other elements 
into account. The two horizontal lines indicate the sizes of the major and minor axes of the synthesized 
beam. 
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Fig. 8. High-resolution (uniform weighting) CO contours at the top of an unsharp-masked HST F814W 
image. The contours are 0.3, 0.4. 0.5, 0.6, 0.7, 0.8, 0.9 times the peak intensity, 0.9 Jy/beam*km/s. 
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Fig. 9. Left: Gas orbits in a bar in the rotating frame with the bar (clockwise rotation). The dashed 
hne indicates the bar direction. Right: An orbit in the non-rotating reference frame. The gas is rotating 

2 

clockwise in both panels. We used the potential ^{R,Lp) = (1 — ecos2(/3)^ log(l + (i?/a)^), where a = 1, 
vq = 1 and e = 0.04. The pattern speed of the bar is fib = 0.1. The two inner Lindblad resonances and the 
corotation radii are i?iLR = 1.1,2.0, and i?cR = 9-9, respectively. Assuming a bar length of 6 kpc (77") 
and rotational velocity of 160 kms""'^ at the radius of 3 kpc for NGC 4303, the units of the model become 
a = 300pc (3".8), -yo = 160kms-\ and f^b = 53kms-i kpc"^ 
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Fig. 10. Left: Gas orbits in a bar projected on the sky. The position angle and inclination are —45° and 
30°, respectively. We first rotated the gas orbits (Figure 9 left) counterclockwise by 4°, so that the bar 
(dashed line) runs vertically, like that in NGC 4303. Middle: Density map, calculated from the gas orbits 
based on the speed of gas motions. The radial density profile is set to be exp(— 7-/4) for r > 2 and constant 
for r < 2. We do not show low density regions. The solid lines indicate the major axis (P. A. = —45°). 
Right: Velocity field, calculated from the gas orbit model. 
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Fig. 11. Position- velocity diagrams of the model in Figure 10 observed from the three different position 
angles, i.e. 0° {top), —45° (middle), and —90° {bottom). The diagram with PA= —45° corresponds to 
NGC 4303. The same gas distribution shows quite different patterns; our determination of PA fits the 
data. The solid curves are the rotation curve given by the assumed potential. 
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Fig. 12. Isochroncs of gas flows. The thick sohd hnes indicate 0, 13, 26, 39, 52 Myr after passage over 
the turning points (the offset ridges). Gas orbits (same as Figure 9) are also shown for a reference. 




Fig. 13. Model color evolution of stellar clusters. 
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Table 1. Parameters of NGC 4303. 



Parameter 




Value 


Reference 


Alias 




M61 




Hubble type 




SAB(rs)bc 


1 


Nuclear activity 




H 


2 


Distance (Mpc) 




16.1 


3 


Linear scale of 1" (pc) 




78.1 




P.A.(isophotal) 








Inclination (isophotal) 








Vsys (kms^i) 




1569 


1 


1)25 (arcmin) 




6.46 


1 


(mag) 




10.12 


1 


•500(1-0) (45") (Jykms- 




494 ±87 


4 


'S'co(i-o)(13") (Jykms^ 




193 ±6 


5 



References: (1) de Vaucouleurs, G. et al. 1991; (2) Ho et al. 1997; (3) Ferrarese et al. 1996; (4) Kcnney, 
Young 1988, 45" (FHWM) beam centered at aigso = 12'n9"^21^4 and 81950 = +4°44'58".0; (5) Nishiyama, 
Nakai 2001, 16" (FHWM) beam centered at ^1950 = 12hl9'"21^4 and 61950 = +4°44'58".0 



Table 2. Observational parameters. 



Parameter Value 

Year 1999 Dec. - 2000 Feb. 

Field center 

"1950, -51950 12hl9"^2P.60, +4°45'03".0 

"2000, '52000 12h2r"54^97, +4°28'24".9 

Field of view 65" 

Array configuration AB, C, D 

Observing frequency 114.65925 GHz 

Bandwidth 512 MHz 
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Table 3. Parameters of CO(l-O) cubes. 



Parameter 


Value 


Value 


Configuration 


AB+C+D 


AB+C+D 


Weigliting 


NA 


UN 


Syntiiesized beam 


2". 8 X 1".9,166° 


1".4 X 1".3,11° 


(kms-i) 


10.4 


31.2 


rms (mjybeam""'^) 


21 


24 


Tf, for IJybcam"^ 


17.3 


50.5 


Tb for IJyarcsec"^ 


91.9 


91.9 



Table 4. Kinematic parameters of the nuclear disk 

Parameter Value 
Dynamical center 

ai950, 51950 12M9"^2P.67, +4°45'03".7 

"2000, -52000 12''21"^54^94, +4°28'25".6 

Vsys (kms-^) 1556.5 

P.A. (deg) -44.45 

i (deg) 29.13 
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